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Abstract 

Using gauge/gravity duality, we analytically calculate properties of a strongly coupled striped 
superconductor, with the charge density wave sourced by a modulated chemical potential, in the 
large modulation wavenumber Q limit. In the absence of a homogeneous term in the chemical 
potential, we show that the critical temperature scales as a negative power of Q for scaling dimen- 
sions A < |, whereas for A > |, there is no phase transition above a certain critical value of Q. 
The condensate is found to scale as a positive power of Q such that the gap is proportional to Q. 
We discuss how these results change if a homogeneous term is added to the chemical potential. 
We compare our analytic results with numerical calculations whenever the latter are available and 
find good agreement. 
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I. INTRODUCTION 



High temperature superconductors, or high Tc superconductors, are materials whose su- 
perconducting transition temperature is higher than 30 K. Even though such materials 
were discovered in the late 1980s, a sufficiently accurate theoretical understanding of such 
materials has been lacking (for a review with a strong emphasis on the theoretical issues, 
see for example, Ref. [T]). The difficulty arises from the fact that these high temperature 
superconductors are strongly coupled [2]- With the lack of a general theoretical framework 
in quantum many-particle physics that deals with fermions at finite density, we face the so- 
called fermion sign problem, in which we cannot rely on brute force techniques using lattice 
models to help us solve this problem at strong coupling. However, recently, by borrowing 
an idea that comes from string theory, i.e., gauge/gravity duality, progress has been made 
in our understanding of these strongly coupled superconductors. 

In gauge/gravity duality, the strongly coupled condensed matter systems are mapped to 
a weakly coupled gravitational theory on a spacetime with negative cosmological constant, 
or the so-called anti de-Sitter (AdS) spacetimes. In particular, the study of high Tc super- 
conductor is mapped to a study of Einstein-Maxwell-scalar theory living in four-dimensional 
AdS black holes and the superconducting phase transition is understood as the (unstable) 
black holes forming scalar 'hairs.' For a review on applications of gauge/gravity duality 
in condensed matter physics, see Ref. [31 and references therein. For analytic studies of 
homogeneous holographic superconductor, see for example, Refs. [U [5] . 

Further complication that arises in the high temperature superconductors, such as 
cuprates and iron pnictides, is the competing orders that are related to the breaking of 
the lattice symmetries. At ffist, these orders seem to be unrelated to superconductivity, 
however, a study of the effect of inhomogeneity of the pairing interaction in a weakly coupled 
BCS system [6] as well as numerical studies of Hubbard models [71 [8] suggest that inhomo- 
geneity might play a role in high T^. superconductivity. Furthermore, the recent discovery 
of transport anomalies in La2_a;Baa;Cu04, which are particularly prominent for x = 1/8 [H], 
might be explained under the assumption that this cuprate is a "striped" superconductor 
|10j . Other studies using mean-field theory have also shown that unlike the homogeneous 
superconductor, the striped superconductor exhibits the existence of a Fermi surface in the 
ordered phase [HI [12] and its complex sensitivity to quenched disorder [10] . 
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Therefore, it is of great value to have comprehensive models that incorporate not only 
strong coupling, but also inhomogeneity. The goal will be to have the competing orders 
emerging dynamically, however, in this article, we will follow Ref. [13] where the inhomo- 
geneity is introduced via a modulated chemical potential, with wavenumber Q. Further work 
on the physical properties of our model is needed to determine its applicability to striped 
superconductors reahzed in Nature. 

We are particularly interested in studying the behavior of the system at large Q. For 
non-vanishing homogeneous part of the chemical potential, the behavior of the leading terms 
of the observables are identical to those of a homogeneous superconductor with some small 
Q-dependent corrections. In this paper, we calculate these subleading terms and we find 
that they are not exponentially suppressed as previously suggested, but instead exhibit a 
power law behavior. When the homogeneous part of the chemical potential vanishes, we 
find a very different qualitative behavior compared to the homogeneous case. In particular, 
the critical temperature vanishes when the anomalous dimension of the condensate is larger 
than 3/2. For values of the dimension below 3/2, we show that the critical temperature 
scales as a negative power of Q whereas the condensate scales as a positive power of Q such 
that the gap is proportional to Q. This indicates that, perhaps, the underlying mechanism 
is not based on a correlation length as with BCS superconductorsj^ It would be interesting 
to further study the properties of these superconductors, including two and higher-point 
correlation functions to better understand the role of the correlation length. 

Our discussion is organized as follows. In Section |Tl} we set up the equations that govern 



the system. In Section III, we calculate the critical temperature in the large Q limit. In 



Section IV, we discuss the behavior of the condensate including subleading terms. Finally, 
in Section IVl we summarize our conclusions. 



II. SET-UP 



The minimal requirement to study strongly coupled superconductor using gauge/gravity 
duality is to have a scalar field and a U{1) gauge field living in a spacetime with negative 
cosmological constant. The scalar field is dual to a scalar order parameter of the super- 
conductor, i.e., the condensate, while the U{1) gauge field is dual to the current in the 

^ We thank Stefanos Papanikolaou for a discussion on this point. 
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condensed matter system. 

To study the strong coupling regime of the superconductor, we only need to study the 
gravity theory at the classical level. In particular, we are interested in finding solutions to 

the classical equations of motion whose boundary values are related to the parameters of 
the superconductor. 

The equation of motion for the scalar field is 

^ Da{./^99'''D,^) + ~^'m) = 0, (1) 



while the equation for the U(l) gauge field is 

^ da (v^F"^) = [^*D,^ - , (2) 



V-9 

and Einstein's equations are 



1 3 

^ab — -j^9abR — ^fl'ab = SttGn 



FacFb'^ — ^9abF'^Fcd 

gabDa^{D,^r + [D,^{D,^r + a O 6] - 9abVm)\ ■ (3) 



Here, Da = da — iqAa, where a is the space-time index. 

In the following, we will consider the potential to be quadratic, with the coefficient m to 
be the mass of the scalar field. Furthermore, we will neglect the backreaction of the scalar 
field and the gauge field onto the metric. This is valid when scale invariant quantities, such 
as Gn/J'^, where /i is the chemical potential, are small. As usual, this can be achieved by 
sending Gat to zero while keeping everything else fixed. 

Therefore, the problem of studying a strongly coupled superconductor with an effective 
two spatial dimensions is reduced to studying the dynamics of a scalar field with mass m, 
coupled to a U{1) gauge field in the background of a S+l-dimensional AdS Schwarzschild 
black hole with the following metric 

ds^ = -r^ h(r) dt^ + —— + dx"^ , (4) 
H h[r) 



with 

in units in which the AdS radius is unity (L = 1). The radius of the horizon is r+ and 
the Hawking temperature, which is equivalent to the temperature of the strongly coupled 



system, is given by 

In this coordinate system, the boundary is at infinity. 

By fixing the gauge, we can set Ax = Ay = = with non-zero electrostatic scalar 



0, (7) 

.// / h' 2\ ^, ( A^ m2\ dl^ 



then { 


j;iven 


by 






dlA 








r\h 






dl^ 




' ^ j 





h z ^ 

where we have set g = 1 and performed a coordinate transformation 

z = ^, (9) 
r 

under which h{r) becomes h = 1 — z^, the boundary is at z = 0, while the horizon is at 
z = 1. Here ' denotes a derivative with respect to z. 

The boundary behaviors of the fields are related to the observables in the strongly coupled 
theory as follows 

^ 0) = - — z, (10) 



^ 0) = ^, (11) 



{Oa) 

where /i is the chemical potential, p is the charge density and {Oa) is the expectation value 
of the condensate of the strongly coupled system. As stated before, r+ is related to the 
temperature via Eq. 

To study the effect of inhomogeneity, we are going to consider the following electrostatic 
potential 

A{z, x)= fij^ ^^"^ ^^"^ (-2) cos nQx, (12) 

n>0 

with 

^5('^) = 1, (13) 

n>0 

subject to the boundary conditions A^"\0) = 1 and A^^\l) = 0. 

This set-up is equivalent to the charge density wave (CDW) being sourced by a modulated 
chemical potential. We are particularly interested in the regime where Q is larger compared 
to the temperature scale. 



III. AT THE CRITICAL TEMPERATURE 



At the critical temperature and above, the order parameter vanishes, and thus \1/ = 0. 
Therefore, 

^<-'"-"^Q:f' =0. (14) 

rih 



For n = 0, the solution is 

= 1-z, (15) 

while for n > 0, we can solve the equation by first rescaling the radial coordinate z = 
and then expanding the equation of motion formally in We obtain 

„3 



a?AH=(l + -^i3 + ...)AH, (16) 



which is an expansion around the boundary z = 0. This expansion is only valid for z ^ 
which covers a small part of the entire range of z if Q is large compared to the temperature. 
However, for z > the potential is exponentially small, so the only physically significant 



range of z is the one in which the formal expansion (16) is valid. 
Indeed, the leading-order solution is 

nO sinh^(l-z) 

= cosh z - coth — sinh z = ^^-^^ . (17) 

r+ sinh =2 ^ ^ 

r+ 

The first-order corrections are of order (nQ/r^)'^ and therefore negligible at large Q for 



^ ^ ^5 whereas for z > the potential (17) clearly vanishes exponentially. As noted in 



Ref. [131 this analytical expression agrees well with the numerical solution even for low Q. 



In fact, as Q — > 0, (17) becomes exact (A^"^ — )■ 1 — z). 

In the following, we shall concentrate on the case with only two Fourier modes. In other 
words, 5^"^ = for n > 1, and therefore, let us set 5^^^ = 5 and 5^°^ = 1 — 5. 

Let us define 

^ = ^^^^ 1_ J2 i^^"^ (z) cos nQx, (18) 
with F*^°)(z = 0) = 1, such that the expectation value of the condensate is 

((^a) = (Oa)^"^ cos nQx, (19) 

n>0 

with 

(Ca)^"^ = (Ca)^°^ F^''\z = 0). (20) 



Here, A = A± = | ± ^/ 1 + is the dimension of the condensate operator Oa- We shall 
examine the range 



2<A<3. 



(21) 



where A > | (A < |) for A = A_|_ (A = A_), corresponding to masses in the range 



> > -|. 



A. Variational Method 



The equation of motion for the scalar field can be written as 



^^2A~2 pin)" ^ (/,^2A-2y pin)' ^ _ ^)^2 _ ^2g2^2] 



,2A-4 pin) 



,2A-2 



(22) 



where /ic (Qc) is the chemical potential (wavenumber) in units of the horizon radius at the 
critical point, 

/xc = ^ , Qc = — . (23) 
Looking at the n = equation, the eigenvalue fic minimizes the expression 

2 _ /o dz z2A-4 {/,^2(^(0)')2 _ (/^ _ ^^^2 (i7^(0))2| 

~ dz ^F(o){(l - 5)2^(0)^(0) - 2(5(1 - 5)AWAWFW + f ^(i)' (FW + F(2))} ' 

(24) 

Since at large Q, it is expected (and numerical study confirms) that F(o) ^ F^^\ F^'^\ the 
above expression reduces to 



dz ^(F(0))2{(1 - 5)2^(0)2 + f A(i)'} 

22A-4 {/,^2^^{0)')2 _ _ 1)^2 (i7^(0))2| 

dz^ii- 5Y Am\Fmy + ^^^^^ 



(25) 



In the absence of a homogeneous term in the chemical potential (5 = 1), Eq. (25) becomes 

22A 



r(2A 



.Qf-i / rf^z2^-^{/i^2(F(°)V-(/^-l)A'(F(°))2} 

~ 1) 



(26) 



Using the trial function F^"^ = 1 — az^, we obtain 



22A /2A2-3A + 6 



r(2A - 1) V 2(2A + 1 
The minimum is attained at 

A2(2A + 1) 



ATI"+2 1^ 



2A-1 

c 



a 



which results in the estimate 



i2A-l 



(A + 1) (2A2-3A + 6)' 
A 2A2 + 9A + 6 



Qi 



2A-1 



(27) 



(28) 



(29) 



r(2A - 1) (A + 1)2(2A2 - 3A + 6) 
Recall that for the above derivation to be valid, we needed Qc ^ 1 (i-e., Q > T^). Therefore, 
/ic > 1. Since 

O O 3-2A 

^ (30) 



,, 2A-1 



it follows that for A > ^, Q/ fi < 1, suggesting that there is a critical value of Q above which 
there is no phase transition. This sudden drop in the critical temperature has been observed 
in the numerical studies, see Fig. [TJ 

T 




0.284 



0.288 // 



FIG. 1. Critical temperature Tc for (5 = 1 and A = 2 obtained numerically. For A > 3/2, there 
is a critical modulation Q* , at which Tc discontinuously drops to zero. Here, Q* = 0.288//. The 
lower branch in Q < Q* regime is unstable. 



For A < |, we deduce the estimate of the critical temperature. 



T 



3 

Air 



i2A-l 



A 



2A2 + 9A + 6 



F(2A - 1) (A + 1)2(2A2 - 3A + 6) 



2A-3 /Q 



(31) 
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For A = 1, we have /i^ = 1.7 Qc and 

Tc _ 0.14 

^ " Q^- ^^^^ 
Thus the critical temperature has a power law behavior in the large Q limit. Of course, this 
is only an estimate of Tc. Nevertheless, the exponent of Q is confirmed by the perturbative 
method to be discussed later, which yields not only the exponent but also the multiplicative 
coefficient accurately in the case 6 = 1. We shall also compare this with the numerical 
results to assess the accuracy of the results of the variational method. 



In the presence of a homogeneous term in the chemical potential {6 ^ 1), Eq. (25 ) becomes 



j^' dz ^ {1 - 5)Hi - z)HF(^)y 

dz ^2^-4 {hz^F^'^yf -{h- 1) A2 r(2A - 1) 5^ 



(1-5)2 (1-^)2 (F(0))2 ' ^c^^ 



(33) 



where we have approximated the second integral in the denominator by evaluating the 
integrand near the horizon. We note that the subleading term is not exponentially suppressed 
as suggested, based on numerical calculations [I3], but is a power law in Q. This is to be 
contrasted with the 1/ logQ behavior seen in the weak coupling BCS calculation [6J. 

Since the second term in the last line is subleading, we can find the eigenvalue fic by first 
minimizing the leading term using a trial function 

= 1 - az\ (34) 

For A = 1, the minimum is attained at a ~ 0.24, which yields 



The critical temperature is 



2 _ 1-27 _ 0.94^2 



^.0.21(1-.) + ^. (36) 
/i Q/fi 



We note that in obtaining the last line of Eq. (33), we have assumed that the first term 
of the denominator is much greater than the second one. Obviously, this assumption is not 
always valid. In particular, it is easy to see that the limit 5 — )■ 1 does not commute with the 
large Q limit. 



B. Perturbative Method 



To gain more insight, let us solve the equation of motion for F*^°^ by treating the electro- 
static potential as a perturbation to a leading order solution, which is nothing but that of a 
scalar field on an AdS Schwarzschild black hole without any background U{1) gauge field 

^.«)% + !(f^) _ ^^<») ^ 0. (37) 

Here, Fq°^ is the leading term of the solution. Strictly speaking, this perturbation is valid 
only when 5 ~ 1, however, as we shall see in a bit, the result is a good approximation even 
when 5 is far away from unity 

The solution that satisfies the correct boundary conditions at 2; = is 



where 2F1 is the Gauss hyper geometric function. There is another solution, which corre- 
sponds to the solution with the correct boundary conditions for A — )■ 3 — A 

^(o)_.3-2A z. /^3-A 3-A_2(3-A)_3 



Using perturbation theory, we obtain the next leading order solution, which is given by 



f p{o)(^\ r dz' ~(o) A (0) p(o)/ X r dz' (0) A 

3-2A^'^° ^'^Jo ^'2-2A^0 ^0 [^)J^ ^,2-2A^0 ^ 







3-2A 

where 



[F^'\z)~a(z)-F^'\z)a{z)), (40) 

^(.) = (1-5)^A(°)^ + ^. (41) 

As both F^^^ and Fq^^ diverge logarithmically at the horizon, we obtain the following 
singularity for the full solution 

F(o)(z ^ 1) = fS'^z ^ 1) + Fi''\z ^ 1) + • • • 



fa log(l — z)' 

where 



2(3-A) 



r(f) 

r2(|) ' 3-2AJ r2(3^) 3-2A 



(42) 



22A g2A-l ■ 
2 
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a{l)^{l-Sra^^^ + l^„ (43) 



with 



dz 



.(0) 



z 



2-2A 



dz 



z 



2-2A 



h 
h 



7(0) 
' 



.(0) 



(44) 



In obtaining Eq. (43), we have approximated the subleading integrals by evaluating the 



integrand near the boundary 2; = 0. 

Requiring regularity at the horizon, we obtain 

p f' 2(3-A) 



1 



3-2A \^r(f) r^(stA) 
(!) 



2(3--A) 
3 



.(0) 



r(f) 12(3^) 




r(2A 



3-2A r(f) r2(3^) 22A(3-2A)g 



2A-1 * 



(45) 



where we have dropped the term proportional to l/Q^ on the last line. Both terms in Eq. 



(45) are small, but to connect this result with the result from variational method Eqs. (35) 



and (31), we make two different approximations in which one is a lot larger than the other. 



When the first term is significantly larger than the second one, for A = 1, we have 

2 1.19 0.9152 

/^c = 7 



5y 



5YQc 



(46) 



which agrees well with Eq. (35). We would like to emphasize that this agreement is valid 



even when 5 is far away from unity. 

When the second term is significantly larger than the first, we have 



r (f ) {^) 22^(3 - 2A) Q 



2A-1 
c 



r2 (f) Y (^ 2(3-A) j r(2A- 1) 

For A = 1, (5 = 1, we then have //^ = 1.55 Qc and 

0.15 



52 



(47) 



T 

^ c 



(48) 



which agrees with the previous result obtained by variational method Eq. (32). 



We would like to note that for 5 = 1 and A > 3/2, Eq. 47 does not have any solutions 



This is related to the fact that for 5 = 1 and A > 3/2, there is a critical Q* such that 
Tc{Q > Q*) = 0, as we have mentioned earlier. 



11 



Let us also comment on the A = 1/2 unitarity limit, which is singular. To approach it, 



we introduce a cutoff A in Q-space. From Eq. (|45|), we have 

p ,'2(3-A)^ 



1 



r(2A - 1^ 



{Q 



1-2A 



A 



1-2A 



(49) 



f,l{Q) /x2(A) r(f) T^'-^) 22A(3-2A) 

If A > 1/2, we can safely take the limit A — )■ oo, in which A^"^"^ — )■ and /Uc(A) has a finite 
limit. This is not so for the case of A = 1/2. Taking the limit A — )■ 1/2, we have 

1 1 



lim 



r(|) n(|)4 ^°^g' 



(50) 



./i2(g) /i2(A)_ 

which shows that the limit A — t- oo is not well defined and that for A = 1/2, the chemical 
potential is no longer a physical quantity. A well defined physical quantity would be 



d 



dlogQ 



r^i) r(i) 6^ 



(51) 



0.10 




T 



5 = 0.3 



10 

Qo 



FIG. 2. The critical temperature as a function of Qc for A = 1. The black lines are obtained by 
numerical analysis and the dashed red lines, which are closer to the numerical results, are obtained 
by perturbative method. 



Lastly, let us compare the critical temperatures obtained by the variational and per- 
turbative methods with the numerical result. This is depicted in Fig. |2} The agreement 
between the results obtained by analytical methods and the numerical results for Qc ^ 3 is 
remarkable. We also show the critical temperature as a function of Q/fi for different values 
of 6 in Fig. [3j 
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FIG. 3. The critical temperature as a function of Q/ ^ for A = 1. The black solid lines, dashed 
red lines and the dashed blue lines are obtained by numerical analysis, perturbative method and 
variational method, respectively. 



C. Higher Modes 

Using the perturbative method, we can also study the behavior of the subleading term 
of the condensate for 1/2 < A < 3/2. For 5 7^ 1, the n = 1 mode does not vanish and at 
sub-leading order, it is given by 



^ 

1 ° [(Oa)(°) 3-2A Jo ^'2-2A-o ^ -0 

+ 3-2A io ° ' ^^^^ 

where Fq^^ and Fq^^ obey 

d. [h ^2^-2 - (zA^ + ^2A-2 ^(1) ^ (53) 

with the boundary conditions Fj^\z = 0) = 1 and F^'' ^ z^'"^^. With /i^ given by Eq. (45), 
the boundary value 

is fully determined by demanding regularity at the horizon. To determine this, let us intro- 
duce the tortoise coordinate, defined as 

.X dz' 
r, = (3 - 2A) 



2A-2 



hz' 

3-2A^/3 3-2A^, 

— - — Biz""- — - — ,0 , (55) 
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where B {x; a, b) is the incomplete Euler beta function. We note that near the boundary 
r^,{z — 7- 0) = z^~'^^. In this tortoise coordinate, the equation of motion for Fq^^ and F^^^ can 
be rewritten in the Schrodinger form 



dr^X'^ + VF^" = 



(1) 



(56) 



with the "potential" V given by 



V 



(3 - 2Ay 



Qc 



4(A-1) 



n^-^^ exp -rt'^ 



3 - 2A, 

Rescaling the tortoise coordinate f* = Ql"^^ r*, at large Qc, we have 

\ 2 4(A-1) 



(1) 



~ 3-2A p(l) 



3 - 2A 

whose solutions are given in terms of modified Bessel functions 

TT 



0, 



7W 




2^ 



r(^^) sin--^ 



and 



Q 



3-2A 



(57) 



(58) 



(59) 



(60) 



These solutions satisfy the boundary conditions FQ^\r^ = 0) = 1 and FQ^\r^ = 0) = f* but 
they diverge at the horizon. Their ratio at large Qc is then given by 



hm ^ = ii^lii^Z 23-2A 



-1 I 



(1) 



sm 







7r 



7r(3-2A) 



Q 



2A-3 
c 



(61) 



Substituting this into Eq. (52) and demanding the full solution to the n = 1 mode be regular 
at the horizon, we obtain 

-, 



Q' 



l + O 



1 



(gA)W_ ^(i-^) r(^)r(^)r(2A-i) . / (3-2A) \ 

(Oa)(°) 3-2A 7r22A-3 2 ; 

(62) 

The leading order coefficient is plotted in Fig. [4j As this vanishes for A = 1 and /ic(A = 1) 
is finite, the leading order is 0{fil/Ql). 
Similarly, we obtain 



(Oa) 



(2) 



52 



(Oa)(°) 4(3 -2A) 



r(5z|A) Y(^) r(2A-l) 



7r22A-i 



sm 



j^^(3-2A)j 1 









1 + 



(63) 
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(Oa)(") \^5{l-5) 



0.2 






0.8 




1.2 


1.4 


-0.2 










-yi 











A 



FIG. 4. The ratio |0^|(o) x^5(i-5) ^ ^ function A. 



The dependence of this ratio on A is depicted in Fig. |5j Again, we see that for A = 1, the 
bracketed quantity vanishes and the ratio is of order 0{iil/Ql). 
It is easy to see that for 5 7^ 1, when A 7^ 1 



(Oa) 



(n) 



/i 



2[n/2] 



where [n] denotes the largest integer < n, and 



(64) 



(Oi) 



(n) 



(Ol)(0) 



o 



2[n/2] rp[n/2] ' 



Q3[n/2 



(65) 



for A = 1. For S = 1, the odd modes vanish while the result for the even modes stays the 
same. This result confirms that at large Q, higher modes are negligible. 

IV. BELOW THE CRITICAL TEMPERATURE 



Below the critical temperature, \E' 7^ and we have to include its contribution toward the 
electrostatic potential. From Eqs. [7] and 18, and by neglecting the higher modes, we have 



A"'> 



n 



rih 



15 



0. 



(66) 



\ 0.03 










0.02 










0.01 










-0.01 


0.8 


1.0\^ 


1.2 


1.4 


-0.02 











FIG. 5. The ratio |^^||o) j?fg ^ ^ function A. 



Rescaling z = ^ z, we have 



in aM^) ?2(a-i) 



where we have dropped terms of order 



A. A = 1 



A 



(67) 



The A = 1 case is a particularly easy one to solve and the solution is given by 

sinh (1 - ,) 



sinh 



(68) 



Substituting this into the equation of motion for the scalar field and demanding regularity 
at the horizon, we have 



rl _ (I) 
/.^ 3r3(|) 



a(l) -a(l) 



(69) 



where a(l) and a(l) are defined in Eq. (40), but now with (n = 0, 1) given by mSn. 
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For 5 — l,we obtain 



(0)2 



Multiplying by /i^, and solving for the condensate, we obtain 



(70) 



{0.r = Q\l{Yj-^, (71) 

showing that in the absence of a homogeneous term in the chemical potential, the condensate 
is proportional to the wavenumber Q. 

Including a homogeneous term {S ^ 1) changes this behavior. The case when the tem- 
perature is not too far away from the critical value, such that the expectation value of the 
condensate is small, will be treated when we study the case of general A, but here let us 
assume that the condensate is large (low temperature regime). Then 



r. r3 (I) 
3r3(i) 



(1 - 5Y 5^ 



(72) 



Multiplying by /i^, we obtain for 5 1 and (Oi)^^^ < Q, 
Therefore, 

{O,)^0)^6A^ + O(l/Q) . (74) 
For Q ^ we similarly obtain at leading order 

(0,)(«)=6.4L_?^^. (75) 



B. Near Tc 

Let us now consider the case of general A. For n — 0, the solution that satisfies the 
boundary conditions is given in terms of the modified Bessel function of the second kind 
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For n > 0, let us first consider a temperature not too far away below T^, such that (Oa) 



(0) 



is small and the right hand side of Eq. (67) can be treated as perturbation. The solution is 
then expanded as 



^(«) ^ ^(n) ^ ^(n) ^ 



(77) 



where 



(n) 



sinh^(l - z) 



(n) 



sinh 



2 



(0) 



(n) 



1 ^ sinh ^z' 

dz' 



A 



(n) 



rf-^nQ I ° Jo {z>y('-^)^' 



1 sinh^(2'-;z 

(^/)2(1-A) 



A, 



(n) 



,2A-1 



-nQz/r+ 



(Oa) 



2(2A - 1) nQ rl^~^ 



(78) 



For small (Oa)^"-*, the zero mode of the electrostatic potential can also be expanded as 



^(0) = l-z- 



(2A - l)22A+i _ + l)z2A + 2z (Ca) 



^2A 



(79) 



(2A-1)2A(2A + 1) 

It is easy to see that by substituting A = 1 into the above equations, we recover the result 



of Eq. (68) in this regime. 



Substituting this into the equation of motion for the scalar and demanding regularity at 



the horizon, we obtain Eq. (69), but now with 



A 



1 sinh^ ^(1 - 2) 

2 sinh' ^ 

r+ 



,2A-1 



-2Qz/r^ 



(Oa) 



2(2A - 1) Q 



„2A-1 ' 



(80) 



for (5 = 1. We deduce 

rl _ (f ) r 



,2A-1 



/i2 r(f) r2(3^) 22A(3 _ 2A) g2A-i 



r(2A- 1) 



r(4A - 2) 1 (Oa)^°^' 

2A - 1 22A-1 g2A 



(81) 



Mulitplying by /i^, after some algebra, we obtain the condensate 
(Oa)^'^ = 



'(2A- l)22A-ir(2A-l) 



r(4A - 2) 

showing that for 5 = 1, in the large Q limit. 



3-2A 



(82) 



(Oa)'/^ocQ. 
18 



(83) 



Similarly, in the case 5 7^ 1, we have 

(!) 



2(3-A) 
3 



+ 



where o'q^ and o^q' are given by Eqs. 

dz 



r(f) r^(^) 
(t) ^ 



2(3-A) 
3 



r(2A - 1) 52 



2(3~A) 
3 



(0) ~(0) 



r (f ) n (^) 



44 



3-2A r(^) r2(3^) 22A(3_2A)g 

(1-^)^ (0a)(°)' 

(2A - 3) A (2A - 1) (2A + 1) rl^ 

and 



2A-1 



a? 



z 
^ dz 



1-2A 



(0) (2A-1)^ 



(2A + l)^2A-i ^ 2 



^2 + 2 + 1 



2-2A 



(2A - 1)22A _ + l)z 

z'^ + Z + 1 



2A-1 



7(0) 
' 



Mulitplying the above by /i^, we get for the condensate 

(Oa)(°) 



7T\ 1 



J^2 



7^2 

c 



Q 



2A-1 



where 



7 



f ) a/A(4A2 - 1) 





(f) 


r 


^2(3-A) 
^ 3 




f) 


r2 


rn 


r2 


(t) 


r 


^2(3-A) 
^ 3 



(0) _ ~{0) 



1/2 



3-A 



a 



(0) ~{o) 



(84) 



(85) 



(86) 



(87) 



The expectation value of the condensate for 5 = 1 and A = 1 obtained analytically is 
compared to the numerical results in Fig. [6j 




FIG. 6. The expectation value of the condensate for 6 = 1 and A = 1 for Qc = 6, 33 and 55. The 
solid black lines and the dashed blue lines depict the numerical and analytical results, respectively. 
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We can also obtain the higher modes of the condensate using a similar calculation as 



the one in Section III C but with the electrostatic potential given in Eqs. 79 and 78 The 
results are 



((9a)(°) 3-2A 
and 



■ _ r(^)r(3^)r(2A-i) 



l + O 



(88) 



(Oa) 



(2) 



6' 



(Ca)(°) 4(3 -2A) 



r(5^) r(3^) r(2A-i) 



7r(3-2A) 



,2--T- - ij q5 [l + O (^ 

(89) 

As before, for A = 1, the (9(/i^/(5^) terms vanish, and the n = 1, 2 modes are suppressed 
by a factor of an order 0(/i^ T/Q^) compared to the n = mode. For 5 = 1, the odd modes 
for the condensate vanish. 



Lastly, let us consider a low temperature where the right hand side of Eq. (67) is a lot 
larger than unity. We then have 



^(,) ^ (0a)W-"2- 



(90) 



r(^)4A^ ■ 

Substituting this into the equation of motion for the scalar and demanding regularity at the 



horizon, we again obtain Eq. (69), but with 



A 



2(5 + 



2 ; P (^) A^ r+ ^ 



Arf 



(91) 



Furthermore, since | < ^ < 1 is not an integer, we can express the modified Bessel function 
in terms of exponential and hypergeometric functions. The integral a(l) then becomes 



a(l) 



1-2(5 + 



1-2(5 + 



2 

352 



(iz 1 



' Arf 



,2-2A 



2A-1 1-A 
2 A A A 



c.(0) ^ M 1.. 1. 2(0a)(°)^^ 
^0 1^1 I 2 2A' A' Ar^ 



2A-1 



(0a)(°) 
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(92) 



where we have approximated the integral by evaluating the integrand at 2; = 0. Similarly, 
we find a(l) to be of order {{Oa)^^^) , which is negligible compared to a(l). Therefore, 



1-25 + f r^(f)r 



2(3-A) 
3 



2A-1 



/i^ 3-2A r(f) r2(3^) 2^^^ V((^a)(°) 

Multiplying by /i^, for 5 = 1, we obtain the condensate 

1-1/A 



(93) 



1 



(0)1 VA 



2 A2A-1 



r(2A - 1) 



2A-1 



3-2A 
2A-1 



Q 



(94) 



showing again that the gap, {O^y^^, is proportional to Q. 

For 5 7^ 1, similarly, we obtain the condensate at low temperature 

1 



(0)1 VA 



7T, 



3-2A 
J" \ 2A-1 



T 



1 + 



Q 



2A-1 



(95) 



where 



Air 



7 



3 (2A^)i/A 



(1-25 + 



r2 


rA^ 


r 


^2(3-A)^ 


< 


2A-1\ 
A ) 




r 






-A 
3 


) 



3-A\ "0 



(0) ~(0) 







(96) 



where a'^^ and d'^^ are again defined in Eqs. (44). 

We do not have a comparative plot for low temperature regime due to the fact that the 
numerical solution to the non-linear field equations become increasingly cumbersome as the 
temperature approaches zero. 

By requiring regularity of the higher modes of the scalar field at the horizon, we can obtain 



the information concerning the higher modes of the condensate, as was done in Section III C 
The results are 

2-A 

2 



(Oa) 



(1) 



and 

(Oa)(^) 
It is easy to see that 



2-A 



4(3 -2A) 2! VA 



//Q \(n) / 2[n/2] 

^^^^ oi [(Oa)(°T~~ 



(Oa)(°) 

where [n] denotes the largest integer < n, but with the odd modes vanishing for 5 = 1 



(97) 



(98) 



(99) 
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V. CONCLUSION 



In this article, we have discussed the properties of a striped holographic superconductor in 
the probe limit at large modulation wavenumber Q. We calculated the critical temperature 
Tc and the expectation value of the condensate (0a) below Tc analytically for arbitrary 
values of the scaling dimension A. 

We found that in the absence of a homogeneous terms in the chemical potential, both 
Tc and (Oa) have a power law behavior for large Q. In particular, the critical temperature 
behaves as 

TcOcQ-i^, (100) 
while the power of the condensate is such that the gap 

{O^y/'^o^Q. (101) 

We also found that the odd modes of the condensate vanish, while the higher even modes 
are suppressed 

^^^''-'s^'fCV (102) 



(c»a)(°) - VQ'"/ 

In the case in which a homogeneous term is included in the chemical potential, both 
Tc and (Oa) approach constant values in the large Q limit, but the subleading terms are 
powers of Q. These constant values are the corresponding values for the homogeneous 
superconductors with chemical potential We also found that the higher modes of the 
condensate are suppressed 

'°^'"'<-0{^^- (103) 



(Oa)(°) ~ VQ^^"/^' / ' 

The behavior of the gap indicates that the underlying mechanism does not rely on a 
correlation length, unlike in the weakly coupled BCS superconductors. The role of the cor- 
relation length is worth exploring further by calculating transport coefficients and correlation 
functions in the large Q limit. Work in this direction is in progress. 

In the weak coupling regime, including fluctuations results in a dip at the small Q regime 
in which critical temperature decreases as Q decreases pass a certain critical value j6j|. We 
expect to see similar behavior once backreaction is included in the holographic calculation. 
Work in this direction is also in progress. 

Lastly, we would also like to note that in the absence of homogeneous term in the chemical 
potential, the zero mode of the condensate is non- vanishing. In other words, the 5 = 1 case 
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does not correspond to the so-called "pure" pair density wave (PDW) state. It will be 
interesting to find a holographic relaization of this pure PDW state and study its properties. 
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